Spiking GABAergic OPC tumor cells prolong survival in IDH1 mutant glioma.

Prior studies have shown that glioma cells form synapses with neurons to receive synaptic inputs. To discern if glioma cells can send outgoing electrochemical signals in the form of action potentials (APs), we employed Patch-sequencing on surgically-resected human glioma slices. Results showed that half of patched cells in IDH1 mutant (IDH1mut) tumors demonstrate select properties of both neurons and glia and fire single, short APs. To define the transcriptional profiles of these hybrid cells (HCs) and discern if they are tumor derived, we developed a computational tool, Single Cell Rule Association Mining (SCRAM), to annotate features in each cell individually. SCRAM revealed that HCs represent a heterogenous group of tumor and non-tumor cells that are uniformly defined by both GABAergic neuron and oligodendrocyte precursor cell (GABA-OPC) transcriptional signatures. We found that GABA-OPC tumor cells express requisite voltage-gated ion channels needed to fire APs. We validated our findings in human single cell and bulk RNA-seq datasets, confirming that GABA-OPCs represent 40% of IDH1mut tumor cells, correlate with survival outcomes in IDH1mut human patients and are also found in select molecular subtypes of IDH1 wild-type tumors. These studies describe a new cell type in IDH1mut glioma with unique electrophysiological and transcriptomic properties.

Targeted gene expression profiling predicts meningioma outcomes and radiotherapy responses.

Surgery is the mainstay of treatment for meningioma, the most common primary intracranial tumor, but improvements in meningioma risk stratification are needed and indications for postoperative radiotherapy are controversial. Here we develop a targeted gene expression biomarker that predicts meningioma outcomes and radiotherapy responses. Using a discovery cohort of 173 meningiomas, we developed a 34-gene expression risk score and performed clinical and analytical validation of this biomarker on independent meningiomas from 12 institutions across 3 continents (N = 1,856), including 103 meningiomas from a prospective clinical trial. The gene expression biomarker improved discrimination of outcomes compared with all other systems tested (N = 9) in the clinical validation cohort for local recurrence (5-year area under the curve (AUC) 0.81) and overall survival (5-year AUC 0.80). The increase in AUC compared with the standard of care, World Health Organization 2021 grade, was 0.11 for local recurrence (95% confidence interval 0.07 to 0.17, P < 0.001). The gene expression biomarker identified meningiomas benefiting from postoperative radiotherapy (hazard ratio 0.54, 95% confidence interval 0.37 to 0.78, P = 0.0001) and suggested postoperative management could be refined for 29.8% of patients. In sum, our results identify a targeted gene expression biomarker that improves discrimination of meningioma outcomes, including prediction of postoperative radiotherapy responses.

COLLAGENE enables privacy-aware federated and collaborative genomic data analysis.

Growing regulatory requirements set barriers around genetic data sharing and collaborations. Moreover, existing privacy-aware paradigms are challenging to deploy in collaborative settings. We present COLLAGENE, a tool base for building secure collaborative genomic data analysis methods. COLLAGENE protects data using shared-key homomorphic encryption and combines encryption with multiparty strategies for efficient privacy-aware collaborative method development. COLLAGENE provides ready-to-run tools for encryption/decryption, matrix processing, and network transfers, which can be immediately integrated into existing pipelines. We demonstrate the usage of COLLAGENE by building a practical federated GWAS protocol for binary phenotypes and a secure meta-analysis protocol. COLLAGENE is available at https://zenodo.org/record/8125935 .

Enhancing Fairness in Disease Prediction by Optimizing Multiple Domain Adversarial Networks.

Predictive models in biomedicine need to ensure equitable and reliable outcomes for the populations they are applied to. Unfortunately, biases in medical predictions can lead to unfair treatment and widening disparities, underscoring the need for effective techniques to address these issues. To enhance fairness, we introduce a framework based on a Multiple Domain Adversarial Neural Network (MDANN), which incorporates multiple adversarial components. In an MDANN, an adversarial module is applied to learn a fair pattern by negative gradients back-propagating across multiple sensitive features (i.e., characteristics of individuals that should not be used to discriminate unfairly between individuals when making predictions or decisions.) We leverage loss functions based on the Area Under the Receiver Operating Characteristic Curve (AUC) to address the class imbalance, promoting equitable classification performance for minority groups (e.g., a subset of the population that is underrepresented or disadvantaged.) Moreover, we utilize pre-trained convolutional autoencoders (CAEs) to extract deep representations of data, aiming to enhance prediction accuracy and fairness. Combining these mechanisms, we alleviate biases and disparities to provide reliable and equitable disease prediction. We empirically demonstrate that the MDANN approach leads to better accuracy and fairness in predicting disease progression using brain imaging data for Alzheimer's Disease and Autism populations than state-of-the-art techniques.

Federated generalized linear mixed models for collaborative genome-wide association studies.

Federated association testing is a powerful approach to conduct large-scale association studies where sites share intermediate statistics through a central server. There are, however, several standing challenges. Confounding factors like population stratification should be carefully modeled across sites. In addition, it is crucial to consider disease etiology using flexible models to prevent biases. Privacy protections for participants pose another significant challenge. Here, we propose distributed Mixed Effects Genome-wide Association study (dMEGA), a method that enables federated generalized linear mixed model-based association testing across multiple sites without explicitly sharing genotype and phenotype data. dMEGA employs a reference projection to correct for population-stratification and utilizes efficient local-gradient updates among sites, incorporating both fixed and random effects. The accuracy and efficiency of dMEGA are demonstrated through simulated and real datasets. dMEGA is publicly available at https://github.com/Li-Wentao/dMEGA.

Single-Cell Profiling of CD8+ T Cells in Acute Myeloid Leukemia Reveals a Continuous Spectrum of Differentiation and Clonal Hyperexpansion.

Comprehensive investigation of CD8+ T cells in acute myeloid leukemia (AML) is essential for developing immunotherapeutic strategies beyond immune checkpoint blockade. Herein, we performed single-cell RNA profiling of CD8+ T cells from 3 healthy bone marrow donors and 23 newly diagnosed (NewlyDx) and 8 relapsed/refractory (RelRef) patients with AML. Cells coexpressing canonical exhaustion markers formed a cluster constituting <1% of all CD8+ T cells. We identified two effector CD8+ T-cell subsets characterized by distinct cytokine and metabolic profiles that were differentially enriched in NewlyDx and RelRef patients. We refined a 25-gene CD8-derived signature correlating with therapy resistance, including genes associated with activation, chemoresistance, and terminal differentiation. Pseudotemporal trajectory analysis supported enrichment of a terminally differentiated state in CD8+ T cells with high CD8-derived signature expression at relapse or refractory disease. Higher expression of the 25-gene CD8 AML signature correlated with poorer outcomes in previously untreated patients with AML, suggesting that the bona fide state of CD8+ T cells and their degree of differentiation are clinically relevant. Immune clonotype tracking revealed more phenotypic transitions in CD8 clonotypes in NewlyDx than in RelRef patients. Furthermore, CD8+ T cells from RelRef patients had a higher degree of clonal hyperexpansion associated with terminal differentiation and higher CD8-derived signature expression. Clonotype-derived antigen prediction revealed that most previously unreported clonotypes were patient-specific, suggesting significant heterogeneity in AML immunogenicity. Thus, immunologic reconstitution in AML is likely to be most successful at earlier disease stages when CD8+ T cells are less differentiated and have greater capacity for clonotype transitions.

Integrated clinical genomic analysis reveals xenobiotic metabolic genes are downregulated in meningiomas of current smokers.

INTRODUCTION: Meningiomas are the most common primary intracranial tumor. Recently, various genetic classification systems for meningioma have been described. We sought to identify clinical drivers of different molecular changes in meningioma. As such, clinical and genomic consequences of smoking in patients with meningiomas remain unexplored. METHODS: 88 tumor samples were analyzed in this study. Whole exome sequencing (WES) was used to assess somatic mutation burden. RNA sequencing data was used to identify differentially expressed genes (DEG) and genes sets (GSEA). RESULTS: Fifty-seven patients had no history of smoking, twenty-two were past smokers, and nine were current smokers. The clinical data showed no major differences in natural history across smoking status. WES revealed absence of AKT1 mutation rate in current or past smokers compared to non-smokers (p = 0.046). Current smokers had increased mutation rate in NOTCH2 compared to past and never smokers (p < 0.05). Mutational signature from current and past smokers showed disrupted DNA mismatch repair (cosine-similarity = 0.759 and 0.783). DEG analysis revealed the xenobiotic metabolic genes UGT2A1 and UGT2A2 were both significantly downregulated in current smokers compared to past (Log2FC = - 3.97, padj = 0.0347 and Log2FC = - 4.18, padj = 0.0304) and never smokers (Log2FC = - 3.86, padj = 0.0235 and Log2FC = - 4.20, padj = 0.0149). GSEA analysis of current smokers showed downregulation of xenobiotic metabolism and enrichment for G2M checkpoint, E2F targets, and mitotic spindle compared to past and never smokers (FDR < 25% each). CONCLUSION: In this study, we conducted a comparative analysis of meningioma patients based on their smoking history, examining both their clinical trajectories and molecular changes. Meningiomas from current smokers were more likely to harbor NOTCH2 mutations, and AKT1 mutations were absent in current or past smokers. Moreover, both current and past smokers exhibited a mutational signature associated with DNA mismatch repair. Meningiomas from current smokers demonstrate downregulation of xenobiotic metabolic enzymes UGT2A1 and UGT2A2, which are downregulated in other smoking related cancers. Furthermore, current smokers exhibited downregulation xenobiotic metabolic gene sets, as well as enrichment in gene sets related to mitotic spindle, E2F targets, and G2M checkpoint, which are hallmark pathways involved in cell division and DNA replication control. In aggregate, our results demonstrate novel alterations in meningioma molecular biology in response to systemic carcinogens.

Single-Cell Profiling of CD8+ T Cells in Acute Myeloid Leukemia Reveals a Continuous Spectrum of Differentiation and Clonal Hyperexpansion.

Comprehensive investigation of CD8+ T cells in acute myeloid leukemia (AML) is essential for developing immunotherapeutic strategies beyond immune checkpoint blockade. Herein, we performed single-cell RNA profiling of CD8+ T cells from 3 healthy bone marrow donors and 23 newly diagnosed (NewlyDx) and 8 relapsed/refractory (RelRef) AML patients. Cells co-expressing canonical exhaustion markers formed a cluster constituting <1% of all CD8+ T cells. We identified two effector CD8+ T cell subsets characterized by distinct cytokine and metabolic profiles that were differentially enriched in NewlyDx and RelRef patients. We refined a 25-gene CD8-derived signature correlating with therapy resistance, including genes associated with activation, chemoresistance, and terminal differentiation. Pseudotemporal trajectory analysis supported enrichment of a terminally differentiated state in CD8+ T cells with high CD8-derived signature expression at relapse or refractory disease. Higher expression of the 25-gene CD8 AML signature correlated with poorer outcomes in previously untreated AML patients, suggesting that the bona fide state of CD8+ T cells and their degree of differentiation are clinically relevant. Immune clonotype tracking revealed more phenotypic transitions in CD8 clonotypes in NewlyDx than in RelRef patients. Furthermore, CD8+ T cells from RelRef patients had a higher degree of clonal hyperexpansion associated with terminal differentiation and higher CD8-derived signature expression. Clonotype-derived antigen prediction revealed that most previously unreported clonotypes were patient-specific, suggesting significant heterogeneity in AML immunogenicity. Thus, immunologic reconstitution in AML is likely to be most successful at earlier disease stages when CD8+ T cells are less differentiated and have greater capacity for clonotype transitions.

Targeted gene expression profiling predicts meningioma outcomes and radiotherapy responses.

Background: Surgery is the mainstay of treatment for meningioma, the most common primary intracranial tumor, but improvements in meningioma risk stratification are needed and current indications for postoperative radiotherapy are controversial. Recent studies have proposed prognostic meningioma classification systems using DNA methylation profiling, copy number variants, DNA sequencing, RNA sequencing, histology, or integrated models based on multiple combined features. Targeted gene expression profiling has generated robust biomarkers integrating multiple molecular features for other cancers, but is understudied for meningiomas. Methods: Targeted gene expression profiling was performed on 173 meningiomas and an optimized gene expression biomarker (34 genes) and risk score (0 to 1) was developed to predict clinical outcomes. Clinical and analytical validation was performed on independent meningiomas from 12 institutions across 3 continents (N = 1856), including 103 meningiomas from a prospective clinical trial. Gene expression biomarker performance was compared to 9 other classification systems. Results: The gene expression biomarker improved discrimination of postoperative meningioma outcomes compared to all other classification systems tested in the independent clinical validation cohort for local recurrence (5-year area under the curve [AUC] 0.81) and overall survival (5-year AUC 0.80). The increase in area under the curve compared to the current standard of care, World Health Organization 2021 grade, was 0.11 for local recurrence (95% confidence interval [CI] 0.07-0.17, P < 0.001). The gene expression biomarker identified meningiomas benefiting from postoperative radiotherapy (hazard ratio 0.54, 95% CI 0.37-0.78, P = 0.0001) and re-classified up to 52.0% meningiomas compared to conventional clinical criteria, suggesting postoperative management could be refined for 29.8% of patients. Conclusions: A targeted gene expression biomarker improves discrimination of meningioma outcomes compared to recent classification systems and predicts postoperative radiotherapy responses.

XCVATR: detection and characterization of variant impact on the Embeddings of single -cell and bulk RNA-sequencing samples.

BACKGROUND: RNA-sequencing has become a standard tool for analyzing gene activity in bulk samples and at the single-cell level. By increasing sample sizes and cell counts, this technique can uncover substantial information about cellular transcriptional states. Beyond quantification of gene expression, RNA-seq can be used for detecting variants, including single nucleotide polymorphisms, small insertions/deletions, and larger variants, such as copy number variants. Notably, joint analysis of variants with cellular transcriptional states may provide insights into the impact of mutations, especially for complex and heterogeneous samples. However, this analysis is often challenging due to a prohibitively high number of variants and cells, which are difficult to summarize and visualize. Further, there is a dearth of methods that assess and summarize the association between detected variants and cellular transcriptional states. RESULTS: Here, we introduce XCVATR (eXpressed Clusters of Variant Alleles in Transcriptome pRofiles), a method that identifies variants and detects local enrichment of expressed variants within embedding of samples and cells in single-cell and bulk RNA-seq datasets. XCVATR visualizes local "clumps" of small and large-scale variants and searches for patterns of association between each variant and cellular states, as described by the coordinates of cell embedding, which can be computed independently using any type of distance metrics, such as principal component analysis or t-distributed stochastic neighbor embedding. Through simulations and analysis of real datasets, we demonstrate that XCVATR can detect enrichment of expressed variants and provide insight into the transcriptional states of cells and samples. We next sequenced 2 new single cell RNA-seq tumor samples and applied XCVATR. XCVATR revealed subtle differences in CNV impact on tumors. CONCLUSIONS: XCVATR is publicly available to download from https://github.com/harmancilab/XCVATR .

Privacy-aware estimation of relatedness in admixed populations.

BACKGROUND: Estimation of genetic relatedness, or kinship, is used occasionally for recreational purposes and in forensic applications. While numerous methods were developed to estimate kinship, they suffer from high computational requirements and often make an untenable assumption of homogeneous population ancestry of the samples. Moreover, genetic privacy is generally overlooked in the usage of kinship estimation methods. There can be ethical concerns about finding unknown familial relationships in third-party databases. Similar ethical concerns may arise while estimating and reporting sensitive population-level statistics such as inbreeding coefficients for the concerns around marginalization and stigmatization. RESULTS: Here, we present SIGFRIED, which makes use of existing reference panels with a projection-based approach that simplifies kinship estimation in the admixed populations. We use simulated and real datasets to demonstrate the accuracy and efficiency of kinship estimation. We present a secure federated kinship estimation framework and implement a secure kinship estimator using homomorphic encryption-based primitives for computing relatedness between samples in two different sites while genotype data are kept confidential. Source code and documentation for our methods can be found at https://doi.org/10.5281/zenodo.7053352. CONCLUSIONS: Analysis of relatedness is fundamentally important for identifying relatives, in association studies, and for estimation of population-level estimates of inbreeding. As the awareness of individual and group genomic privacy is growing, privacy-preserving methods for the estimation of relatedness are needed. Presented methods alleviate the ethical and privacy concerns in the analysis of relatedness in admixed, historically isolated and underrepresented populations. SHORT ABSTRACT: Genetic relatedness is a central quantity used for finding relatives in databases, correcting biases in genome wide association studies and for estimating population-level statistics. Methods for estimating genetic relatedness have high computational requirements, and occasionally do not consider individuals from admixed ancestries. Furthermore, the ethical concerns around using genetic data and calculating relatedness are not considered. We present a projection-based approach that can efficiently and accurately estimate kinship. We implement our method using encryption-based techniques that provide provable security guarantees to protect genetic data while kinship statistics are computed among multiple sites.

SVAT: Secure outsourcing of variant annotation and genotype aggregation.

BACKGROUND: Sequencing of thousands of samples provides genetic variants with allele frequencies spanning a very large spectrum and gives invaluable insight into genetic determinants of diseases. Protecting the genetic privacy of participants is challenging as only a few rare variants can easily re-identify an individual among millions. In certain cases, there are policy barriers against sharing genetic data from indigenous populations and stigmatizing conditions. RESULTS: We present SVAT, a method for secure outsourcing of variant annotation and aggregation, which are two basic steps in variant interpretation and detection of causal variants. SVAT uses homomorphic encryption to encrypt the data at the client-side. The data always stays encrypted while it is stored, in-transit, and most importantly while it is analyzed. SVAT makes use of a vectorized data representation to convert annotation and aggregation into efficient vectorized operations in a single framework. Also, SVAT utilizes a secure re-encryption approach so that multiple disparate genotype datasets can be combined for federated aggregation and secure computation of allele frequencies on the aggregated dataset. CONCLUSIONS: Overall, SVAT provides a secure, flexible, and practical framework for privacy-aware outsourcing of annotation, filtering, and aggregation of genetic variants. SVAT is publicly available for download from https://github.com/harmancilab/SVAT .

The evolving privacy and security concerns for genomic data analysis and sharing as observed from the iDASH competition.

Concerns regarding inappropriate leakage of sensitive personal information as well as unauthorized data use are increasing with the growth of genomic data repositories. Therefore, privacy and security of genomic data have become increasingly important and need to be studied. With many proposed protection techniques, their applicability in support of biomedical research should be well understood. For this purpose, we have organized a community effort in the past 8 years through the integrating data for analysis, anonymization and sharing consortium to address this practical challenge. In this article, we summarize our experience from these competitions, report lessons learned from the events in 2020/2021 as examples, and discuss potential future research directions in this emerging field.

Evaluation of vicinity-based hidden Markov models for genotype imputation.

BACKGROUND: The decreasing cost of DNA sequencing has led to a great increase in our knowledge about genetic variation. While population-scale projects bring important insight into genotype-phenotype relationships, the cost of performing whole-genome sequencing on large samples is still prohibitive. In-silico genotype imputation coupled with genotyping-by-arrays is a cost-effective and accurate alternative for genotyping of common and uncommon variants. Imputation methods compare the genotypes of the typed variants with the large population-specific reference panels and estimate the genotypes of untyped variants by making use of the linkage disequilibrium patterns. Most accurate imputation methods are based on the Li-Stephens hidden Markov model, HMM, that treats the sequence of each chromosome as a mosaic of the haplotypes from the reference panel. RESULTS: Here we assess the accuracy of vicinity-based HMMs, where each untyped variant is imputed using the typed variants in a small window around itself (as small as 1 centimorgan). Locality-based imputation is used recently by machine learning-based genotype imputation approaches. We assess how the parameters of the vicinity-based HMMs impact the imputation accuracy in a comprehensive set of benchmarks and show that vicinity-based HMMs can accurately impute common and uncommon variants. CONCLUSIONS: Our results indicate that locality-based imputation models can be effectively used for genotype imputation. The parameter settings that we identified can be used in future methods and vicinity-based HMMs can be used for re-structuring and parallelizing new imputation methods. The source code for the vicinity-based HMM implementations is publicly available at https://github.com/harmancilab/LoHaMMer .

Cross-Talk Between Histone Methyltransferases and Demethylases Regulate REST Transcription During Neurogenesis.

The RE1 Silencing Transcription Factor (REST) is a major regulator of neurogenesis and brain development. Medulloblastoma (MB) is a pediatric brain cancer characterized by a blockade of neuronal specification. REST gene expression is aberrantly elevated in a subset of MBs that are driven by constitutive activation of sonic hedgehog (SHH) signaling in cerebellar granular progenitor cells (CGNPs), the cells of origin of this subgroup of tumors. To understand its transcriptional deregulation in MBs, we first studied control of Rest gene expression during neuronal differentiation of normal mouse CGNPs. Higher Rest expression was observed in proliferating CGNPs compared to differentiating neurons. Interestingly, two Rest isoforms were expressed in CGNPs, of which only one showed a significant reduction in expression during neurogenesis. In proliferating CGNPs, higher MLL4 and KDM7A activities opposed by the repressive polycomb repressive complex 2 (PRC2) and the G9A/G9A-like protein (GLP) complex function allowed Rest homeostasis. During differentiation, reduction in MLL4 enrichment on chromatin, in conjunction with an increase in PRC2/G9A/GLP/KDM7A activities promoted a decline in Rest expression. These findings suggest a lineage-context specific paradoxical role for KDM7A in the regulation of Rest expression in CGNPs. In human SHH-MBs (SHH-α and SHH-ß) where elevated REST gene expression is associated with poor prognosis, up- or downregulation of KDM7A caused a significant worsening in patient survival. Our studies are the first to implicate KDM7A in REST regulation and in MB biology.

Multiple approaches converge on three biological subtypes of meningioma and extract new insights from published studies.

One-fifth of meningiomas classified as benign by World Health Organization (WHO) histopathological grading will behave malignantly. To better diagnose these tumors, several groups turned to DNA methylation, whereas we combined RNA-sequencing (RNA-seq) and cytogenetics. Both approaches were more accurate than histopathology in identifying aggressive tumors, but whether they revealed similar tumor types was unclear. We therefore performed unbiased DNA methylation, RNA-seq, and cytogenetic profiling on 110 primary meningiomas WHO grade I and II). Each technique distinguished the same three groups (two benign and one malignant) as our previous molecular classification; integrating these methods into one classifier further improved accuracy. Computational modeling revealed strong correlations between transcription and cytogenetic changes, particularly loss of chromosome 1p, in malignant tumors. Applying our classifier to data from previous studies also resolved certain anomalies entailed by grouping tumors by WHO grade. Accurate classification will therefore elucidate meningioma biology as well as improve diagnosis and prognosis.

Racial and Socioeconomic Disparities in Patients With Meningioma: A Retrospective Cohort Study.

BACKGROUND: Meningiomas are the most common intracranial neoplasms. Although genomic analysis has helped elucidate differences in survival, there is evidence that racial disparities may influence outcomes. African Americans have a higher incidence of meningiomas and poorer survival outcomes. The etiology of these disparities remains unclear, but may include a combination of pathophysiology and other factors. OBJECTIVE: To determine factors that contribute to different clinical outcomes in racial populations. METHODS: We retrospectively reviewed 305 patients who underwent resection for meningiomas at a single tertiary care facility. We used descriptive statistics and univariate, multivariable, and Kaplan-Meier analyses to study clinical, radiographical, and histopathological differences. RESULTS: Minority patients were more likely to present through the emergency department than an outpatient clinic (P < .0001). They were more likely to present with more advanced clinical symptoms with lower Karnofsky Performance scores, more frequently had peritumoral edema (P = .0031), and experienced longer postoperative stays in the hospital (P = .0053), and African-American patients had higher hospitalization costs (P = .046) and were more likely to be publicly insured. Extent of resection was an independent predictor of recurrence freedom (P = .039). Presentation in clinic setting trended toward an association with recurrence-free survival (P = .055). We observed no significant difference in gross total resection rates, postoperative recurrence, or recurrence-free survival. CONCLUSION: Minority patients are more likely to present with severe symptoms, require longer perioperative hospitalization, and generate higher hospitalization costs. This may be due to socioeconomic factors that affect access to health care. Targeting barriers to access, especially to subspecialty care, may facilitate more appropriate and timely diagnosis, thereby improving patient care and outcomes.

Ultrafast homomorphic encryption models enable secure outsourcing of genotype imputation.

Genotype imputation is a fundamental step in genomic data analysis, where missing variant genotypes are predicted using the existing genotypes of nearby "tag" variants. Although researchers can outsource genotype imputation, privacy concerns may prohibit genetic data sharing with an untrusted imputation service. Here, we developed secure genotype imputation using efficient homomorphic encryption (HE) techniques. In HE-based methods, the genotype data are secure while it is in transit, at rest, and in analysis. It can only be decrypted by the owner. We compared secure imputation with three state-of-the-art non-secure methods and found that HE-based methods provide genetic data security with comparable accuracy for common variants. HE-based methods have time and memory requirements that are comparable or lower than those for the non-secure methods. Our results provide evidence that HE-based methods can practically perform resource-intensive computations for high-throughput genetic data analysis. The source code is freely available for download at https://github.com/K-miran/secure-imputation.